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Abstract 

Nonperturbative picture of strong interacting quark-gluon plasma 
is given based on the systematic Field Correlator Method. Equa- 
tion of state, phase transition in density-temperature plane is derived 
and compared to lattice data as well as subsequent thermodynamical 
quantities of QGP. 



1 Introduction 

The perturbative exploring of quark-gluon plasma (QGP) has some difficul- 
ties in describing the physics of QGP and phase transitions. However, it 
was realized 30 years ago that nonperturbative (np) vacuum fields are strong 
(PP) and later it was predicted ([2]) and confirmed on the lattice ([3]) that 
the magnetic part of gluon condensate does not decrease at T > T c and even 
grows as T 4 at large T ([!]). 

Therefore it is natural to apply the np approach, the Field Correlator 
Method (FCM) ([3]) to the problem of QGP and phase transitions, which 
was done in a series of papers ([6]- [10]). As a result one obtains np equation 
of state (EoS) of QGP and the full picture of phase transition, including an 
unbiased prediction for the critical temperature T c (/i) for different number 
of flavors rif. 
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2 Nonperturbative EoS of QGP 



We split the gluonic field A M into the background field and the (valence 
gluon) quantum field a^. A^ = B^ + a^ both satisfying the periodic boundary 
conditions. 

The partition function averaged both in perturbative and np fields is 

Z(V,T) = (Z(B + a)) B>a (1) 

Exploring free energy F(T,fi) = —T\ii(Z(B))b that contains perturba- 
tive and np interactions of quarks and gluons (which also includes creation 
and dissociation of bound states) we follow so-called Single Line Approxima- 
tion (SLA). Namely, we assume that quark-gluon system for T > T c stays 
gauge invariant, as it was for T < T c , and neglect all perturbative interac- 
tions in the first approximation. Nevertheless in SLA already exist a strong 
interaction of gluons (and quarks) with np vacuum fields. This interaction 
consists of colorelectric (CE) and colormagnetic (CM) parts. The CE part 
in deconfinement phase creates np self-energy contribution for every quark 
and gluon embedded in corresponding Polyakov line. An important point 
is that Polyakov line is computed from the gauge invariant qq (gg) Wilson 
loop, which for np Df interaction splits into individual quark (gluon) con- 
tributions. As for CM part - its consideration is beyond the SLA, because 
as has been recently shown in paper ([11]) strong CM fields are responsible 
for creation of bound states of white combinations of quarks and gluons. 

To proceed with FCM we apply the nonabelian Stokes theorem and the 
Gaussian approximation to compute the Polyakov line in terms of np field 
correlators 
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As a result the Polyakov loop can be expressed in terms of "potentials" 
V\ and Vjj 
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In what follows we use the Polyakov line fit (j9j [TO] ) 
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The free energy F(T) of quarks and gluons in SLA can be expressed as a 
sum over all Matsubara winding numbers n with coefficients L n j und and L^ d - 
for quarks and gluons respectively. For nonzero chemical potential /i one can 
keep Lf und ^ ad j independent of /i, treating np fields as strong and unchanged 
by /i in the first approximation. 

The final formulas for pressure of qgp are ([HI [10] ) 
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where z/ = m q /T and 
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The energy density is e = T 2 -^ {Jf) v an d the speed of sound in plasma 



is c; 



dp 

de 



In Fig. 5 c s is shown calculated with the use of (jHJ), (TTOTl and 



compared to lattice data for /i = from ([12]). No lattice calculations has 
yet been done for sound speed at nonzero baryon density, though our theory 
allows to do that and as is shown in ( [13] ) the result for /x > does not differ 
much from the case fi = 0. 
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Figure 1: Fit © of 
Polyakov line for rif = 
and n/ = 2) (black curves) 
to the lattice data Q12J). 

3 Phase transition 
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Figure 2: Analytic ©, flE) 
and lattice (|12j) curves for 
pressure of QGP with Uf — 
0,2 + 1,3 from ([10J). 



To obtain the curve of phase transition one needs to define pressure Pi in 
the confined phase and Pjj in the deconfined phase, taking into account 
that vacuum energy density contributes to the free energy, and hence to the 
pressure: 

Pi — \ £ vac\ + Phardon, Pi I = \ £ tac\ + (Pq + PglfP^- (H) 

Having formulas for pressure (which contain parameter of Lf un( i(x) ((7j)) we 
may write for the phase transition curve T c (/i): 

/ ( ii-jn,)AG, y" 

UW \ 32(p,+ P!l ) ) {U] 

here Ae vac = \e vac — = (11 — |nj)/32AG 2 - in particular, for the ex- 
pected value of AG2/G2 (stand) « 0.4 one obtains T c = 0.27 GeV (n/ = 
0), 0.19 GeV (rif = 2), 0.17 GeV (n/ = 3) in good agreement with lattice 
data. 



4 Summary 

The EoS of QGP is written, where the only np input is the Polyakov line. 
It should be stressed, that only the modulus of the Polyakov line enters in 
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Figure 3: Analytic and lat- 
tice ([12]) curves for energy 
density of QGP with rif = 
2 + 1 and n f = 3 from ([TDj). 



Figure 4: The phase tran- 
sition curve T c (/i) (in GeV) 
from (J9j) as function of 
quark chemical potential \x 
(in GeV) for = 2 (up- 
per curve) and rif = 3 
(lower curve) and AG 2 = 
0.0034 GeV 4 from ([9]). 
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Figure 5: Sound speed for fi = and rif = 3 (blue dashed curve) compared 
to lattice data from ([12]). 



EoS due to gauge invariance. The phase transition curve T c (//) and speed 
of sound c^(T) are obtained and agree well with lattice data. An important 
point of the work is that the only parameter used to receive the final physical 
quantities from the initial QCD Lagrangian is the Polyakov line taken from 
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lattice data, and is in agreement with analytic estimate for T = T c ([Sj)- 
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